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THEORETICAL ANTISYMMETRIC SPAN LOADING FOR WINGS OF ARBITRARY PLAN FORM 

AT SUBSONIC SPEEDS 1 

By John DeYoung 


SUMMARY 

A simplified lifting-surface theory that includes effects of 
compressibility and spanwise variation of section lift-curve slope 
is used to provide charts with which antisymmetric loading due 
to arbitrary antisymmetric angle of attack can be found jor 
wings haring symmetric plan jorms with a constant spanwise 
sweep angle of the Quarter-chord line . Consideration is given 
to the flexible wing in roll. Aerodynamic characteristics due to 
rolling , deflected ailerons , and sideslip oj wings with dihedral 
are considered . Solutions are presented jor straight-tapered 
wings for a range of swept plan jorms. 

INTRODUCTION 

Reference 1 has been for many years the standard reference 
for estimating the stability and control characteristics of 
wings. The lifting-line theory on which this work was 
based gave generally satisfactory results for straight wings 
having the aspect ratios considered ; however, the use of wing 
sweep combined with low aspect ratio lias made an extension 
of this work desirable. Lifting-line theory cannot ade- 
quately account for the increased induction effects due to 
sweep and low aspect ratio; consequently, it lias been found 
necessary to turn to the more complex lifting-surface theories. 

Of the many possible procedures, a simplified lifting- 
surface theory proposed by Weissinger and further developed 
and extended in reference 2 has been found especially suited 
to the rapid computation of characteristics of wings of 
arbitrary plan form. Comparisons with experiment have 
generally verified the theoretical predictions. In reference 
2, this method lias been used to compute for plain, unflapped 
wings, the aerodynamic characteristics dependent on sym- 
metric loading. The same simplified lifting-surface theory 
can be extended to predict the span loading resulting from 
antisymmetric 2 distribution of the wing angle of attack. 
From such loadings the damping moment due to rolling, 
the rolling moment due to deflected ailerons, and the rolling 
moment due to dihedral angle with the wing in sideslip can 
be determined. A recent publication (reference 3) makes 
use of the simplified lifting-surface theory to find span- 
loading characteristics of straight- tapered swept wings in roll 
and loading due to dihedral angle with the wing in sideslip. 
Experimental checks of the theory for the damping-in-roll 
coefficient and rolling moment due to sideslip w r ere very 
favorable. The range of plan forms considered in reference 


3 is somewhat limited and aileron effectiveness was not 
included. The loading due to aileron deflection normally 
involves excessive labor when computed by means of the 
simplified lifting-surface theory; however, development of 
the theory, presented in reference 4, that deals with flap and 
aileron effectiveness for low-aspect-ratio wings provides a 
means by which the simplified lifting-surface method can be 
used to obtain spanwise loading due to aileron deflection. 

It is the purpose of the present analysis to provide simple 
methods of finding antisymmetric loading and the associated 
aerodynamic coefficients and derivatives for wings with sym- 
metric plan forms limited only by a straight quarter-chord 
line over the semispan. Means will be presented for finding 
quickly the aerodynamic coefficients of span loading due to 
rolling, of span loading due to deflected ailerons, and of 
span loading due to sideslip of wings with dihedral. Flexible 
wings, when the flexure depends principally on span loading 
as in loading due to rolling, can be included in the analysis. 
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NOTATION 

■ /& 2 \ 
aspect ratio I ) 

wing span measured perpendicular to the plane of 
symmetry, feet 
wing chord, feet 3 
aileron chord, feet 3 

mean wing chord feet 3 

/local lift" 


local lift coefficient (J" qc~^) 

i ii + /induced drag\ 

induced drag coefficient ^ (jS ) 

. /rolling moment\ 
rolling-moment coefficient — J 

rolling moment due to rolling ‘Dalian 

rolling moment due to aileron deflection ^ d&y 
per radian 

spanwise loading coefficient for unit rolling moment 

(^) 

scale factor 


1 Supersedes XAC A TN 2140, “Theoretical Antisymmetric Span Loading for W r ings of Arbitrary Plan Form at Subsonic Speeds” by John peYoung, 19M. 

2 The word “antisymmetric” is understood to indicate that a distribution ofloading or angle of attack is equal in absolute magnitude on each half of the 'King 

3 Measured parallel to the plane of symmetry. 


but of opposite sign. 
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factors of loading interpolation function 
spanwise loading coefficient or dimensionless circu- 


lation 


( 2 ‘0 or (bv) 


per radian 

spanwise loading coefficient due to aileron deflection 
P el radian 

wing geometry, compressibility, and section lift- 

curve-slope parameter (^) (c~7/3 

integration factors for spanwise loading due to ailer- 
ons 

Mach number 

arbitrary number of span stations defined by 

77 7T 

V = cos — y— 
m + 1 

rate of rolling, radians per second 
wing-tip helix angle, radians 

coefficient depending on wing geometry and indi- 
cating the influence of antisymmetric loading at 
span station n on the downwash angle at span 
station v 

free-stream dynamic pressure, pounds per square 
foot 

wing area, square feet 

ratio of aileron chord to wing chord 3 

free-stream velocity, feet per second 
induced velocity, normal to the lifting surface, 
positive for downwash, feet per second 
lateral coordinate measured from the wing root 
perpendicular to the plane of symmetry, feet 
section angle of attack at span station v , radians 3 
angle of antisymmetric twist of the elastic wing 
produced by the loading due to rolling, radians 3 
rate of change of wing-section angle of attack with 
control -surface angle for constant section lift 

coefficient 3 

compressibility parameter (\|1 M 2 |) 

angle of sideslip, radians 

dihedral angle measured perpendicular to the plane 
of symmetry, radians 

spanwise circulation, feet squared per second 
angle of deflection of full wing-chord control surface, 
radians 3 

angle of deflection of full- wing-chord control surface, 
measured perpendicular to the hinge line, radians 

dimensionless lateral coordinate (^^2) 


/ G \ 


(%) 

\pbj2Vj’ 

6 

trigonometric 


v /aileron span\ 

dimensionless aileron span ^ 6/2 — / 

spanwise center of pressure on one wing panel 


edge of the aileron span, radians 
ratio of section lift-curve slope at a span station 
2 7T 

v to — / both at the same Mach number 

sweep angle of the wing quarter-chord line, positive 
for sweepback, degrees 
compressibility sweep-angle parameter 


j^tan 1 degrees 

/ tip chord \ 

taper ratio ( , T 1 

\root chord/ 


dimensionless lateral coordinate 


<t> trigonometric spanwise coordinate (cos -1 77), radians 

SUBSCRIPTS 

7 ij v integers pertaining to specific span stations given by 

n tt vir 

Tj — cos - — or 7] — cos c - 

o o 

k pertaining to span station k 

c. p. center of pressure 

a aileron 

t pertaining to fraction-of-wing-chord ailerons 

T wing tip 

lx wing root 

av average or mean 

DEVELOPMENT OF METHOD 

The simplified lifting-surface method used herein replaces 
a lifting surface by a lifting vortex located at the wing one- 
quarter-chord line. The boundary condition for determining 
the vortex strength distribution specifies that, along the 
three-quarter-chord line of the wing, there shall be no flow 
through the lifting surface. In effect, this specifies that, at 
the three-quarter-chord line, the ratio of the velocity nornnd 
to the mean camber line (induced by the bound and trailing 
vortices) to the velocity of the free stream shall equal the 
sine of the angle of attack. 

Span loadings are theoretically additive. Since the sym- 
metric angle -of-attack distribution contributes only to sym- 
metric loading, it follows that the antisymmetric loading is 
independent of sym metrically distributed wing twist or cam- 
ber; hence, to find antisymmetric loading, it is only necessary 
to consider the loading resulting form the antisymmetric 
distribution of the angle of attack across the wing span. In 
the subject case, such a distribution is experienced by the 
wing as induced angle due to rolling/ the effective twist due 
to aileron deflection, or sideslip of the wing with dihedral 


5 Measured parallel to the plane of symmetry. 

* In considering the case of the angle induced by rolling as equivalent loan antisymmetric distribution of twist, it must be noted that account should be taken of the fact that a rolling wing 
leaves a twisted vortex trail; whereas a twisted w ing does not. The difference in induction effects on the wing of the straight and twisted vortex is considered insignificant here, as has been 
assumed in other analyses. 
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For an antisymmetric angle-of-attack distribution, the 
loading distribution will be equal in absolute magnitu e^on 
each semispan, but of opposite sign. The loading therefore 
needs only to be found over the semispan, and, since the 
loading is zero at the wing root, only span stations outboarc 
need be considered. The mathematical development of the 
simplified lifting-surface method for the case of antisym- 
metric loading is given in appendix A. As shown m 
\ (m-lV2 linear equations in terms of loading distribution 
are obtained which satisfy the wing angle-of-attack condi- 
tions* at the three-quarter-chord line at m stations n, wild 
m is an arbitrary' odd integer. These equations are repre- 
sented by the summations 


a,= S?»X, "=1,2,3, 


m- 


CD 


whore 


antisymmetric angle of attack at wing station " 
coefficients that for a given value of m depend on 
wing geometry, compressibility, and section htt- 
cut'vc slope 

loading coefficients at span stations n 
The application in appendix A of the present report is with 
7/1 = 7 Since the loading at the midspan station is hno*n to 
be zero, consideration is required of only three stations: 

n = l 2 3 equal to wing semispan positions of if cos (mri 

= 0.924; 0.707; and 0.383. Equation (1) thus becomes 


shown in appendix B, p n can be conveniently presented as a 
function of two parameters, namely, a compress, hle-sweep- 
angle parameter defined as A^tan" 1 (tan A//3) and a param- 
eter TL involving the ratio of wing span to wing chord and 
variable section lift-curve slope, defined by 

( 3 ) 


"-''■(£) Gjs) 


.'ratio of experimental section lift-curve slope at span 
station v to the theoretical value of 2ir/0, both at t ,e 
same Mach number 
c y wing chord at span station v 
The value <h is a scale factor given by 
d,=0.061 for v=\\ 

= 0.234 for v=2) ( 4 ) 

= 0.381 for r=3 

Equation (3) can be written in alternative form that gives 
77„ in terms of wing geometry parameters that are more 
significant ; thus 

(5) 


O n 


a„ = PvnOn, 
71 — 1 


1.2,3 


( 2 ) 


where the integer v pertains to span station rj cos (vir'S) 
To obtain the loading coefficients f7.= (c i r/26)„, it remains 
only to evaluate the coefficients p, n and the spanvisc 
variation of the antisymmetric angle of attack 

EVALUATION OF COEFFICIENTS p yn 

Since m is chosen, becomes a function only of wmg 
geometry, compressibility, and section lift-curve s ope. 11 
effects of compressibility and section lift-curve slope are 
equivalent to a change in wing plan form 6 and can be ac- 
counted for by a proper adjustment of the p rn values. * . 


whore 


to 27 r/j3 


CvjCa 




TJ V d v ^ [(./ K ar ) (c ^/e ar )^l 

ratio of average section lift-curve slope 
both at the same Mach number 
span wise distribution of section lift-<ut\c slop* 
for a given Mach number 
spanwisc distribution of the wing chord 
compressible aspect ratio and average section 
lift-eurve-slope parameter 

* of equation (5) gives an 


effective 


The distribution of K v jK ar ma^ 


T1, ‘‘ term ( k ,/ k ..)(*/0 

aerodynamic taper of a ^ing. — - - 

vary with Macli number, particularly at transonic spec, s 
(e. i due to spanwisc variation of airfoil section) . However, 
since' the distribution contributes to taper effect, the loading 
distribution and not the total loading will be appreciably 

With TL determined from equations (3) or (5) and (4), the 
values of *»„, nine in all, are presented in figure 1 where p vn is 
given as a function of TL for various values of A*. 

: . . Thi. substitution of a r for sin lias the effect of incroaslng the 

5 The reader should note that the boundary condition was fixed 'Tt^n be^en reTdffy « this 

value of loading on the wmg above that necossap " .. arr assumed to move downst ream in a horizontal plane rom t continue to be satisfied, the strength of the bound 

“““ 

"niscusslon” and Sii the appftvlis n. 
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(0 i—2, n~3, *=0.231. 
Figure 1. — Continued. 
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(ti) „=3, h = 2, 6=0.381; 
(i) *=3, n = 3, 6=0.381. 
Fjgure 1.— Concluded. 
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For tlie case of straight-tapered wings with arbitrary 
section lift-curve-slope distribution for which the chord 
distribution is specified by taper ratio, evaluation of equation 


(5) is given in figure 2 where for each of the three 

(p A I K ac ) 

span stations is shown as a function of taper ratio. 


EVALUATION OF ANTISYMMETRIC ANGLE-OF-ATTACK DISTRIBUTION a , 

The antisymmetric angle- of-a ttack distributions most 
commonly encountered are those resulting from rolling 
wings, aileron deflection, and sideslip of wings with dihedral. 
Evaluation of the angle-of-attack distributions for these 
various cases is outlined in the sections immediately following. 

Rolling wings. — For the case of the rigid wing, the induced 
velocity normal to the wing surface is equal to the up wash 
velocity experienced by the rolling wing. Thus, at span 
station v 



where pb/2V is the tip helix angle. It should be noted that 
the relation given by equation (6) assumes the wing structure 
to be rigid in that the distribution of is completely defined 
by the linear distribution of helix angle. In the case of 
flexible wings, however, the expression for a v must be modified 
to account for the streamwise angle-of-attack change which 
^^piary occur due to bending or torsional deflections. In this 
case, 


OLy=— ??„ + Aa„ 


(7) 


where Aa v represents the modifying influence of flexibility. 
Normally, Aa v is not considered for straight wings since only 
the effect of torsion (which is usually small) is involved. 
On swept wings, however, the effect of bending can cause 
Aa v to be quite large so that the a v distribution may be 
affected considerably. Due to the interaction existing 
between the aerodynamic and structural forces, Aa, cannot 
be determined directly, but must be found through equations 
of equilibrium or by iteration. With the loading for the 
rigid wing provided, however, the iteration procedure 
becomes relatively easy to apply. The first approximation 
of a p is found from the loading of the rigid wing and further 
refinements of a y may be found utilizing the successive 
loadings for the flexible wing as determined. 

Deflected ailerons. — Where the spanwise distribution of 
the angle a v is to be considered equivalent to antisymmetric 
aileron deflection, it must suffer a discontinuity at the span- 
wise end of the control surface. The loading when such a 
discontinuity is present can be duplicated by a proper 
distribution of antisymmetric twist. In appendix C, the 
antisymmetric twist distribution required by the present 
theory to give accurate span loading distribution due to 
ailerons is found with the aid of zero-aspect-ratio wing theory 
given by reference 4. To minimize the computation involved, 
it is convenient to consider both the case of outboard and 
inboard ailerons. 

1. Outboard ailerons . — With m — 7, three different aileron 
spans can be conveniently defined for the outboard ailerons. 
For the aileron spans rja, measured from the wing tip inboard, 
the antisymmetric twist distribution required per unit 
deflection of full-wing-chord ailerons, a,/5, is given by 


Cast 1 

I 

TT 

III 

V* 

0. 169 

0. 444 

0. 805 

a i 

5 

1.003 

0.971 

0.998 

a 2 

a 

.017 

996 

.991 

5 

.006 

.014 

.978 


(8) 


Inboard ailerons. — With m — 7, three different aileron spans 
can be conveniently defined for the inboard ailerons. For 
the aileron spans rja measured from the wing midspan out- 
board, the antisymmetric twist distribution required per unit 
deflection of full -wing-chord ailerons, a„/5 , is given by 


Case 

IV 

V 

VI 

Vo 

0. 550 

0.831 

1.000 

6 

0. 044 

0.013 

1.016 

a i 

$ 

-.017 

.961 

.979 

«3 

6 

I. 087 

1.09.5 

1. 101 


(0) 


Sideslip of wings with dihedral. For calculating the roll- 
ing moment caused by dihedral angle for the sideslipping 
wing, the effect of the skewness of the vortex field in altering 
the effects of the dihedral angle will be assumed to be small 
(as assumed in reference 3). The problem then simplifies to 
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finding the rolling moment due to antisymmetric angle of 
attack with the unskcwcd vortex field. The solution to this 
problem is the same as for the ailerons which has already been 
solved. 

The antisymmetric distribution of angle of attack for the 
sideslipping wing with dihedral is given by 

( 10 ) 

where 

a v effective angle-of-attack distribution 
j3 angle of sideslip measured positive in the counterclock- 
wise direction from the plane of symmetry 
T dihedral angle 

The wing parameter T is not affected by compressibility. 
Equation (H)) is approximate for small values of (3 and l\ 

For unit £r over the span of the ailerons considered, 

pr=5 (ii) 

can be substituted for 5 in equations (8) and (9). 

APPLICATION OF METHOD 

For the cases of antisymmetric angle-of-attack distribu- 
tions resulting from rolling, aileron deflection, or sideslip with 
dihedral, it is possible to present a set of simultaneous equa- 
tions which are required for the solution of the load distribu- 
tion for an arbitrary plan form. With the loading known, 
integration formulas can be given to determine aerodynamic 
coefficients. 

The loading-distribution coefficient G n determined from 
the solutions of the simultaneous equations, are functions of 
p> n which has been shown in a preceding section to be a 
function of wing geometry, compressibility, and section lift- 
curve slope. The aerodynamic coefficients are integrations 
of the load distribution and, therefore, will also be a function 
of wing geometry, compressibility, and section lift-curve- 
slope parameters. Application of the method to the general 
solution for arbitrary chord distribution is outlined and solu- 
tions are presented for the case of straight taper. 

GENERAL SOLUTION 


and 


— 0 . 924 + ^ PnG v + P 12 G 2 + PviG-i 
— 0.707 + ^“-yy=p-2iG l + p 22 G 2 + i) a G 3 ► 
— 0. 3 8 3 + pfioy = Ci + fta ft + ft 3 ft 


( 13 ) 


where G n 


G n 


- 1 and A a, is the incremental angle of attack 
pbj 2v * 


due to aeroelastic effects. 

2. Loading distribution .- — The loading-distribution coef- 
ficient is given by G— Cic/2b. Other forms of the loading 
coefficient are given by the identities 


1 CiC_ Ci CiC 

2 A c av 2 A CiCav 


(14) 


The loading is known to be zero at rj= 0 and 1 and is deter- 
mined at three intermediate span stations. Values of loading 
at other span stations can be obtained from a loading function 
derived in appendix B or, with equations (B23) or (B24) of 
appendix B, the loading can be found at span positions 
i|=0.981, 0.831, 0.556, and 0.195. 

3. Rolling moment . — The da in ping-in-roll derivative for the 
solutions of equations (12) or (13) is derived in appendix B 
and given by 

—=Tr [^+0-707 (ft+ft)] (15) 

Kar 1 0 \ K av / 


4. Span wise center o f pressure . — The equation giving center 
of pressure on the wing semispan is shown in appendix B to be 


_ j SCifKar 

r,c p ■ pA (O.163ft + O.2480 2 +0. 430ft) 

K av 


1 


0.082 



+ 0.124 



+ 0.215 



(16) 


Aerodynamic characteristics due to rolling. — The solutions 
for the aerodynamic effects due to the rolling wing will be 
found and loading, rolling moment, spanwise center of pres- 
sure, and induced drag will be obtained. 

1. Simultaneous loading equations . — The p vn values are 
obtained from figure 1 and table 1 7 with values of II v given by 
equations (3) or (5). 

The simultaneous equations (2), for the rigid and flexible 
wing, respectively, become: 


where 


0.924 — p\\G PxzGz^r ptfGi 
— 0.707=p 2 i&i+P22^+P236 f 3 ► 
0.383 = p3\GA~ pz 2 G 2 -\- PmGz ^ 


G n - 


G n 


pb/ 2 V 


( 12 ) 


5. Induced drag . — The induced drag is derived in appendix 
B and given by 

—^=1 (;f) [ft 2 + ft 2 + ft 2 -^ 0,(0i + ft)] (17) 

Aerodynamic characteristics due to aileron deflection. — 

The solutions for the aerodynamic effects due to ailerons will 
be found for three different spans of outboard and inboard 
ailerons. Cross plots of these data provide curves for 
arbitrary aileron spans. 

1. Simultaneous loading equations . — The p vn values are 
obtained from figure 1 and table I with values of II v given 
by equations (3) or (5). 

(a) Deflected outboard ailerons . — The aileron spans meas- 
ured from the wing tip inboard are given by q a . The simul- 
taneous solution for antisymmetric spanwise loading due to 


7 Values of p v » beyond the scope of figure I are Included in table I. For values of IT, larger than those included in figure 1 and table I, the p» n curves can be obtained from equation (B&) 
which gives the linear asymptotes of the function. • 
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deflection of any of the three following aileron spans can be 
obtained from the appropriate set of the following equations: 


=PnG\ -\-V\2G1 -i-pxzGz 
= p'2\C\ '\~P'>2p2~\~P'Xi^ ( 18 ) 

= : /> 3 \G\ F P^^iF pZzGz 

w here G n =Gjd . 

(b) Deflected Inboard ailerons . — The aileron spans meas- 
ured from the wing midspan outboard are given by 17 *. The 
simultaneous solution for antisymmetric spanwisc loading 
due to deflection of any of the three following aileron spans 
can be obtained from the appropriate set of the following 
equations: 


Case 

I 

II 

III 


0. 169 

0. 444 

0. 805 

6 

1.003 

0. 971 

0. 998 

CTj 

.017 

.996 

.991 


.006 

.014 

.978 

$ 






= p\\G x + PiJiv + pv/j* 

= P 21 G i~\~ P2zG 2 ~r p2^s (1 fl) 
— pA 1 + Pi2p 2 "f pM@3 


where $„=G n /d. 

2. Loading distribution . — The spanwise loading distribu- 
tions due to various aileron configurations include: 

(a) Full -wing-chord ailerons . — The loading is known to be 
zero at 77=0 and 1, and is determined at three intermediate 
span stations. With equation (C13) and tables C6, Bl, and 
C7, the loading can be found at span stations tj= 0.981, 0.831, 
0.556, and 0.195 for each of the aileron spans considered. 
With these given points and the knowledge that the slope of 
the loading distribution curve is theoretically infinite at the 
point of angle-of-attack discontinuity (aileron spanwise end), 
the loading distribution can be faired. 

( b ) Constant fraction of wing-chord ailerons . — The spanwise 
loading of constant fraction of wing-chord ailerons is equal to 
the product of the loading due to full- wing-chord ailerons and 
the effective change of angle of attack with aileron angle, 8 
dajdd. The factor dajdd is a function of the ratio of aileron 
chord to wing chord t=cjc. The change of section angle of 
attack with aileron angle dajdd is presented in figure 3, which 
is reproduced from figure 18 of reference 5. 

Although figure 3 taken from reference 5 limits the Mach 
number range to Mach numbers less than 0 . 2 , this limitation 
is believed to be unwarranted since theory indicates that 
dajdd is unaffected by compressibility for the two-dimen- 
sional wing. However, as indicated in reference 4, da/d5 is 
strongly affected by low’ aspect ratio and will change appreci- 
ably if the parameter pA becomes much less than two; hence, 
the values of dajdd from figure 3 appear to be valid for /Li> 2. 

g In using dajds here, it should be noted that the assumption is made that the effect! v 
airfoil section is taken as being parallel to the plane of symmetry and that the section 
^pproQches a two-dimensional section. The validity of this assumption can he Questioned 
however, limited checks with experiment show it to be at least approximately correct. 


Case 

IV 

V 

VI 


0. 556 

0.831 

I. 000 

« 

0. 044 

0. 013 

1. 016 | 

t 



6 

-.017 

.961 

.979 

6 

1.087 

1.09, 5 

1. 101 


FiornE 3. Variat ion of lift-effectiveness parameter with aileron chord ratio, f—y • Average 
trailing-edge angle about 10°; Af< 0.2. Curves from reference 5. 

(c) Arbitrary spanwise distribution of aileron chord. -Tlu* 
aileron can be divided into several spans with constant 
dajdd , then the total loading is the sum of the products of 
the full-wing-chord loading of eacli span and its respective 

dajdd. 

3. Foiling moment. — The rolling moment can he found for 
the following aileron configurations: 

(a) Full-wi ng-chord ailerons.— The spamvise loading due to 
aileron deflection cannot be integrated with sufficient 
accuracy with equation (15). In appendix C, a similar 
integration formula is developed that applies to each given 
aileron span. Equation (CIO) and table C5 give 

— *=( Qh^i + hM+h^) (20) 

Kac \ Kat / 


where for each of the eases of equations (18) and (19) the 
h n values are given by 


Case 

I 

II 

III 

IV 

V 

VI 

h [ 

0. 140 

0. 139 

0. 138 

0.146 

0. 141 

0. 140 

A] 

.199 

. 196 

. 196 

.200 

. 197 

. 198 

hi 

. 145 

. 139 

.138 

.140 

. 139 

. 140 


(6) Constant fraction of wing-chord ailerons . — For constant 
fraction of wing-chord ailerons with aileron angle measured 
parallel to the plane of symmetry, the aileron effectiveness is 
given by 

cZtt/flCA (21 ) 

K<jp d d \ K av ) 

(c) Arbitrary spanwise distribution of aileron chord. The 
deflection of ailerons for which t varies spanwise on the 
wing can be considered as an equivalent wing-twist distribu- 
tion. The effective antisymmetric twist of the wing is 
given by 

«-=s s ' (22) 
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whore dajdd is now a function of spanwisc position. The 
antisymmetric angle-of-atlack distribution given by equation 
(22) can be divided into spanwisc steps of constant angle of 
attack and the total rolling moment can be found by the 
summation of tbe rolling moment due to each spanwisc step. 
The rolling moments of tbe spanwisc steps are obtained from 
a curve of rolling-moment coefficient (3CijK av as a. function 
of unit antisymmetric angle of attack from tbe wing root 
outboard. This step method is tbe procedure used in 

reference 1. . 

A curve of fiC ijs a v as a function of unit antisymmetric 
angle of attack from the wing root outboard can be obtained 
from the solutions of equation (19) for the cases IV, V, and 
VT. An additional point can be obtained from the solution 
of case III of equation (18), applying the relations (dis- 
cussed later) existing between inboard and outboard ailerons. 
The rolling moment due to the twist, given by equation (22) 
can be obtained, by a method other than the step method, 
from the integral given by 




Kav 


-f 


(]a ^l J K av) 

Is <h 


(23) 


which can be integrated numerically by taking the graphical 
slopes of pC, s lK ttV which is a function of extent of unit anti- 
symmetric angle of attack from tbe wing root outboard. 

' 4. Spanwise center of pressure and induced drag. — Span- 
wise center of pressure and induced-drag integration formulas 
for loading due to ailerons are not given; however, equations 
(16) and (17) can give approximate integrations of the load- 
ing to obtain center of pressure and induced drag. 

5. Additional considerations: 

(a) Relation between aerodynamic characteristics for out- 
board and inboard ailerons— The spanwise loading distribu- 
tions due to outboard amt inboard ailerons bear a simple 
relation to each other. Since loading is linearly propor- 
tional to angle of attack, loadings are directly additive. 
Then, for outboard and inboard ailerons with the spanwise 
ends of the ailerons at tbe same span station, 


G i 


'inboara : 


c 

1 ^inboard 


= 6f(ij =i) — ^outboard 


fy (j 

l &{y a = 1) ^outboard 


^ “inboard ^ ^ “outboard 


(24) 


These relations do not apply for nc.r. and C D . since these 
characteristics are not linearly proportional to loading. 

(b) Differential aileron angles . — The effect of a differential 
between aileron angles can lie taken into account by con- 
sidering the Ci s of each wing panel as one-half the antisym 
metric results of equations (20), (21), or (23). The total 
wing rolling moment is then the sum of the products of CiJ 2 
given by equations (20), (21), or (23) and the angle of de- 
flection of each aileron. Although the total rolling moment 
can be found by this procedure, the spanwisc loading distri- 
bution can be found only approximately by the products of 
the antisymmetric unit loading G/d and the deflection of each 
aileron. However, the loading distribution so found will be 

985590-52 3 


quite accurate since this procedure neglects only the small 
change due to the induced effects of the differentially differ- 
ent opposite wing panels. 

(c) Aileron angles measured perpendicular to the hinge 
line— The relationship between aileron angle measured per- 
pendicular to the aileron hinge line and that measured 
parallel to the plane of symmetry is given by 


tan 5 = 


tan 5 


(25) 


COS A; 

where 

A, sweep angle of the aileron hinge line 
8 angle measured perpendicular to the binge line. 

For constant fraction of wing-chord ailerons on straight- 
tapered wings, A, is given by 


tan A, = lan A 


cj\~ 


4(0.75 — f) 

A 


(m) 


(26) 


where t is the fraction of wing-chord aileron measured from 
the wing trailing edge. 

Aerodynamic characteristics due to sideslip of wings with 
dihedral. The total antisymmetric loading due to sideslip 
can be considered as the sum of that due to dihedral angle 
and that due to zero dihedral angle. For the unswept wing, 
the rolling moment due to sideslip for zero dihedral angle is 
generally considered negligible ; however, for the swept wing, 
this effect can he appreciable. In tbe present report, only 
that part due to dihedral angle will be considered for the 
swept and nonswept wings. 

1. Simultaneous loading equations. The p, n values are 
obtained from figure 1 and table 1 with values of given by 

equations (3) or (5). . 

The simultaneous equations resulting from the substitu- 
tion of 8-]sr (sec equation (11)) and G=GfdV in equations 
(18) and (19) are applicable in the determination of the 
effects of unit outboard or inboard dihedral angle over the 
span of the ailerons considered. 

2. Rolling moment . — Tbe rolling moment due to various 
dihedral angle distributions include: 

(a) Constant spanwise dihedral angle . — For dihedral angle 
constant for the entire wing semispan, the loading is__given 
by the solution of case VI in equation (19) for G=G/pT and 
the rolling moment from equation (20) becomes 

(o. 1 40^ + 0. 198^rk 0 ■ 1 4 OG'a) (27) 

Kav 

(, b ) Gulled wing . — For the gulled wing, solutions of equa- 
tion (19) for 7j=Glpr gives the loading, and the rolling 
moment from equation (20) becomes 


130 P A (h l G l +h:G t +h 3 Gf) 


(28) 




A plot of the results of cases IV, V, and VT gives the extent 
of unit dihedral angle from the wing root outboard. Then, 
for a gulled wing, the total rolling moment equals the sum 
of products of dihedral angle of each span section and the 
rolling-moment contribution of the respective span sections, 
(c) Variable spanwise dihedral angle ..— ff T varies span- 
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vise, the rolling moment can be obtained by integration as 
in equation (23). The. integral becomes 


where 


d{HC r jK av T) 


is the slope of the curve described in 


(Itj 

part (6) above. 

SOLUTION FOR STRAIGHT-TAPERED WINGS 

Charts of aerodynamic characteristics for straight-ta- 
pered wings can be presented in terms of geometric, com- 
pressibility, and average section lift-curve-slope parameters. 
These charts provide a ready means of obtaining data 
directly. 

Aerodynamic characteristics due to rolling. -The applica- 
tion of equation (12) for a constant value of section lift- 
curve slope 0 provides the span wise loadings at span stations 
0.383, 0./07, and 0.924 which are presented in figure 4 for a 
wide range of plan forms. The interpolation formula of 
equation (B24) wall give values of loading due to rolling at 
span stations other than those presented. With equation 
(15), the damping-in -roll coefficients &C 1p /k uv can be obtained 
and are presented in figure 5 for a wide range of plan forms. 

Aerodynamic characteristics due to aileron deflection. - 
The application of equation (19), case III of equation (18), 
and equation (20) provide aileron effectiveness in the coeffi- 
cient form fiCijKav for several aileron spans. In figure 0, 
fiCiJtiav is plotted against extent of unit antisymmetric angle 
of attack from the wing semispan root outboard for a range 
of wing parameters. 

As presented, figure 6 gives directly the effectiveness of 
full-wing-chord inboard ailerons for aileron spans measured 
from the plane of symmetry outboard. The effectiveness 
of full -wing-cliord outboard ailerons for aileron spans meas- 
ured from the wing tip inboard is given by figure 6 directly 
by the relations of equation (24). For full- wing-chord 
ailerons located arbitrarily on the wing semispan, the aileron 
effectiveness can be obtained directly from figure G as indi- 
cated in the following example sketch. 



With the full-wing-chord values given above, the effec- 
tiveness of constant fraction of wing-chord ailerons or ailerons 
of arbitrary spnnwise chord distribution can be found through 
use of equat ions (21) or (23) with thet/cr/tft values of figure 3. 

9 Throughout the figures, *•„ , is the constant span wise section lift curve slope or the average 
of a small variation. For large span wise variations of « that follow the function given by 
equation (BU) developed in appendix B, the parameters j9.1/* a , and X can be replaced by the. 

BA*T 

parameters ( <JB ^-, r x)/(i4-x) an<l X ’ n * s P° ctfvol y* For lar 6 c spanwise variations of * 

that do not follow the curve of equation (BU), the simultaneous equations for the general 
solution can be solved for arbitrary distributions of a. The IT, values can be obtained from 
figure 2. 



Figure 4. 


(a) i)= 0.3827. 

-Variation of loading due to rolling coefficient 
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parameter Ap, degrees, for straight -tapered wings. 
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(b) ij= 0.7071. 
Figure 4. -Continued. 


Aerodynamic characteristics due to sideslip of wings with 
dihedral. The application of equation (19), case ITT of 
equation (IS), but with 5=PT, and G=~G/pr, and the use of 
equation (28) provides rolling moments due to dihedral angle 
for the wing in sideslip. These rolling moments are given 
in the coefficient form QC r J^T which is the same function 
of rj as BC and is presented with in figure 0. 

Figure G with equation (29) will provide the rolling moment 









N 3 


3 - 

6 

5 

J>- 

5 - 

4 \ 


l 2\ 


-M=o 


/3A--0 


TT 


-50 -40 -30 - 20 -10 


| A = i.5 


XIEJ 


0 10 20 30 40 50 60 70 80 


(c) tj = 0.9239. 

Figure 4.~ Concluded. 

due sideslip for any symmetric spanwisc distribution of 
dihedral angle. 

For dihedral angle constant spanwisc, the rolling moment 
is given by the value at 77 — 1 in figure 0. These values for 
constant spanwisc dihedral angle are presented in figure 
7 as a function of aspect ratio for various values of sweep 
angle and taper ratio. 

DISCUSSION 

Effects of plan-form parameters on aerodynamic charac- 
teristics for straight-tapered wings arc shown by plots against 
the various parameters. Compressibility is discussed and 
formulas given for a range of plan forms at sonic speeds. 
Theoretical considerations and experimental comparisons 
indicate the order of reliability of the present theoretical 
results. 

STRAIGHT-TAPERED WINGS 

The spanwisc loading distribution due to rolling for several 
plan forms is presented in figure 8. These curves are the 
result of applying figure 4 and the loading interpolation 
formula of appendix B. The loading coefficient is given as 

Vc -^ — ( fy lC ) make the total loading on the semispan 

{Vc.pJA = 0 \U iCav / 

constant and thus show more clearly the changes of dis- 
tribution due to sweep and taper ratio. Figure S shows 
large changes in loading distribution for the zero tapered 
wing. The effects of sweep are generally as expected, 
namely, that sweepback shifts the loading outboard. 

Effects of plan form on the rolling moment due to rolling is 
shown from cross plots of figure 5 which arc presented in figures 
9 and 10. For higher aspect ratio, figures 4, 9, and 1 0 show the 
marked lowering of rolling moment due to sweep. Figure 



(a) x=o. 
(b) X=0.25. 


Figure 5.— Variation of dam ping-in -roll parameter — with compressible swoop parameter 

K av 

Ap, degrees, for straight-tapered wings. 

9 indicates that for low aspect ratio, the rolling moment 
becomes essentially independent of sweep and taper. The 
taper effects on rolling moment as seen in figure 10 are 
small except for values of taper ratio less than 0.25. 

Typical spanwisc loading distributions due to full-wing- 
cliord aileron deflection are shown in figure 1 1 . These curves 
were faired with the aid of the loading interpolation function 
of appendix C and, at the aileron spanwisc end, care was 
taken to make the slope large. 

Wing geometry effects on aileron effectiveness for full- 
chord outboard partial-span ailerons (with aileron angle 
measured parallel to the plane of symmetry) are given in 
figure 12. The geometry effects on fiCijKa «. arc similar to 
those on the damping-in-roll coefficient. Comparison of 
figure 12 (a) with figure 9 shows that C h approaches the 
zero -aspect-ratio value in the same manner as docs Ci v , 
Figure 13 gives comparative effectiveness of inboard and 
outboard ailerons for swept wings. As sweep increases, the 
difference of effectiveness between inboard and out- 
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(c) X=0.50. 

(d) \ = 1.0. 

Fkjure 5.- -Continued. 


board ailerons decreases showing that inboard ailerons for 
highly swept-baek wings approach the effectiveness of out- 
board ailerons. Since dajdb becomes large rapidly at small 
values of t (fig.3), then, for a given aileron area, narrow full- 
span ailerons for swep t-baek wings. may be more desirable 
than larger-chord outboard ailerons. The relative effects of 
figures 12 and 13 apply equally well for constant fraction of 
chord ailerons, since the data would differ only by a constant 
factor da. ld5. 

COMPRESSIBILITY 

From the three-dimensional linearizcd-compressiblc-flow 
equation, it can he shown that the effects of compressibility 
will be properly taken into account if the longitudinal com- 



ic) X s * 1.5. 

FmrRE 5.- Concluded. 



/3Cf- 

sideslip with dihe drul — ~,i per radian squared, for extent of unit antisymmetric angle or 
attack from the wing root outboard. 
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pononts of a wing plan form are increased by the factor ]/fi. 
Or, alternatively, if the linearized compressible flow equation 
be divided through by 0 2 , then the lateral and vertical com- 
ponents of a plan form are decreased by the factor 0. In 
both cases, the incompressible local lift is increased by the 
factor 1//3 and the compressible local lift coefficient can be 
written as the parameter 

With these relations known, an incompressible theory 
can be made into a compressible theory subject to the 
limitations of the linearized compressible flow equation. 
The geometric parameters of a wing are simply sub- 

sti luted by fiA, A,= tan-' and fib. With local lift 

coefficient given by /3c l7 the dimensionless loading 

becomes The wing-chord distribution remains 

2 $b 2b 

unaltered. 

The sonic speed results of reference 4 can be used as a 
limit point in the present theory for a curve of the variation 
of antisymmetric aerodynamic characteristics with Mach 
number. The following equations apply at the speed of 
sound to plan forms with all points of the trailing edge at 
or behind the upstream line of maximum wing span: 

n 

C h-^2“ 



(C) x=i.o. 

FinrRE 6.— Concluded. 

For outboard ailerons, 

4 

C x =■—- sin 3 0. where 7] a = 1 —cos <9 
6 G 

For inboard ailerons, 

4 

(7^==^ (1 —sin 3 6 ), where 77„=cos 6 

Reference 4 shows that aileron effectiveness at the speed of 
sound is independent of the chord wise location of the ai- 
leron binge line, provided the hinge line remains ahead of 
all points of the trailing edge. 

ACCURACY OF THE SEVEN-POINT SOLUTION FOR AILERONS 

The prediction of aileron effectiveness for given aileron 
spans with wing twist determined by zero-aspect-ratio 
theory at only seven span points to satisfy the boundary 
conditions has been theoretically shown to be sufficient by 
comparing results with the computation of a typical 3.5 
aspect ratio, 45° swept wing with 15 span points satisfying 
the boundary conditions. The process of finding aileron 
spans for the 15-point method was the same as that in 
appendix C. The curves showing the variation of C is with 
aileron span for the 7- and 15-point computations were 

identical. .. 

The solution for the angle-of-attack distribution that 
includes a discontinuity can be compared with the solution 
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for the continuous angle-of-attack distribution by consid- 
ering an aileron such that the angle-of-attack distribution 
is equivalent to that of the rolling wing. The damping- 
in-roll coefficient then can be found by use of equation (23) 
which reduces to the form 


C t = 



dC,, 

dr] 


di] 



and integrating by parts 


CV 



c. 


= 1 


This relation states that C lp is equal to the area between a 
curve of figure 6 and the line of Ci h for rj^l. The curves of 
figure G were found by the simplified lifting-surface theory 
with antisymmetric twist determined by zero-aspect-ratio 
theory. The values of C ip obtained in this manner from 
figure G v r ere identical to the C tp values given by simplified 
lifting-surface theory for continuous linear antisynrmetric- 
twist distribution. 

As further theoretical check, the values of rolling moment 
due to constant spamvise dihedral angle are obtained from 
15-point computations in reference 3 for taper ratio equal to 
one, with which the present theory for the 7-point method 
is in exact agreement. 




(a) A p=60°. 

(b) A 0 = 45°. 

(c) A 0 = 0°. 

(d) A(j= —45°. 

Figure 8 . -Span wise loading due to rolling of wings with various taper-ratio and sweep-angle 

BA 

parameters of aspect-ratio parameter — =4. The curve for A =0 serves as a basis for com- 
parison and the factor j c t p • gives the curves constant area equal to 3.454. 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 

Tlie electro-magnetic analogy method of reference 6 
provides damping-in-roll coefficients for an aspect-ratio 
range of unswept, tapered wings. The results of the present 
theory and those of reference 8 are compared in figure 14. 
Except for the taper ratio effects on C Xp the comparison is 
good. The rounded-wing- tip values of C ip given by NACA 
Rep. 635 (reference 1) are included in figure 14. Since 
rounded wing tips generally give values of C* p about 6 
percent lower than straight wing tips, the values of NACA 
Rep. 635 appear to be appreciably too high for lower-aspect- 
ratio wings. The present theory and the theory of refer- 



Fioure 9. Variation of the damping-in-roll parameter - with aspect-ratio parameter 
— f or various sweep angles and taper ratios. 



fiCi 

Fiot're 10— Variation of the damping-in -roll parametc with ta|>or ratio for wings 

with various aspect-ratio and sweep-angle parameters. 


ence 6 approach the value given by the zero-aspect-ratio 
theory of reference 4 quite satisfactorily. The results 
of the present theory may be further assessed by the com- 
parison with the results of low-speed experiment as given 
in figure 15 for the range of plan forms presented. For 
further experimental verification of the accuracy with which 
Ci can be determined by the present theory, the reader is 
referred to reference 3 which supports the theory as well 
or better than figure 15 of the present report. 




Fmt’RK 12.- -Aileron rollmg-moment-coefTicient parameter — — 7 per radian, for outboard 

ailerons as a function of aspcct-rat-io parameter, taper ratio, and eompressiblc-sweep param- 
eter. (a) * = 0.5. 


R 1 

trouRE 13. Aileron rolling-moment parameter — — 9 per radian, for ~ = 4, \=0.5,asafuno 


tion of outboard and inboard aileron spans. 
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Fuii'KE 14.— Comparison of damping-in-roll coefficient C lp of the present theory with those 
of the theories of references 1, 4, and 6. 

The loading distributions due to rolling as given by the 
present theory are compared in figure 16 with low-speed 
experimental results for a range of swept wings. The sweep 
angle seems to have considerably more influence on loading 
distribution as given by experiment than the theory indicates. 
The experimental pressure data, however, were very erratic 
and no firm conclusion can be made. 

Experimental values of rolling effectiveness due to aileron 
deflection are compared with theoretically predicted values 
in a correlation diagram given by figure 17. Included are 
the results of a wide range of plan forms which do not vary 
consistently with an}" geometric parameter or aileron con- 
figuration. Sketches of the plan forms and ailerons are 
drawn about the points of correlation. The theory makes 
use of the curve of figure 3 giving dajdb for a sealed -gap 
aileron over a range of deflection of ±10°. Experimental 
results for aileron deflections greater than 15° measured per- 
pendicular to the hinge line were not included. The correla- 
tion points of figure 17 scatter appreciably; however, the 
mean line of the points does approximate the line of perfect 
correlation. 

Figure 17 does not account for effective plan-form change 
due to dajdb. Only the effectiveness of the low-aspect-ratio 
triangular wing of figure 17 is exceedingly in error, which is 
the result of neglecting plan-form change. 

The plan-form change due to dajdb can, in part, be con- 
sidered analogous to that due to section lift-curve-slope 
change. Thus, the total section lift of a wing, the chord of 
which is reduced by dajdb and which is at an angle of attack 
5, is equal to the lift of the wing-aileron section for which 
the aileron only is deflected at the angle b . This change in 
plan form, unlike the section lift-curve-slope change for 


which the chordwise loading remains constant, does not ac- 
count for a large change in chordwise loading. If the lifting 
line is considered to be at the chordwise center of pressure, 
then, for partial-wing-span ailerons, the lifting line is in 
effect broken at the aileron spamvise end and the present 
theory becomes invalid. For the case of full-wing-span ai- 
lerons, the lifting line in effect remains unbroken and lies 
along the center of chordwise pressure. For this case the 
wdng chord can be reduced by dajdb to account for plan-form 
change; however, although in the limit of zero aspect ratio 
the results are the same as those of reference 4, this procedure 
does not with sufficient accuracy account for the chordwise 
loading shifting aft at intermediate aspect ratios. For con- 
trol surfaces, the effective plan-form change due to dajdb is 
appreciable for the low-aspect-ratio wings such that in the 
limit of zero aspect ratio the spamvise loading is independent 
of the ratio of aileron chord to wing chord (reference 4). 
However, for moderate aspect ratios, dajdb can be used with- 
out accounting for plan-form changes as comparison with 
experiment indicates. 

Experimental values of C^/r are not compared with the 
present theory since reference 3 gives ample support of the 
theory. 

CONCLUDING REMARKS 

The determination of antisymmetric loading for arbitrary 
wings is shown to be easily obtained by the solution of three 
simultaneous equations. The coefficients of the simulta- 
neous equations are presented in charts of parameters that 
include wing geometry, compressibility, and section lift- 
curve slope as arbitrary quantities. Thus the loading for 
an arbitrary antisymmetric angle-of-attack distribution can 
be simply found once the angle-of-attack distribution is 
chosen. 

For the important cases of antisymmetric loading, roll, 
and aileron deflection, the angle-of-attack distribution is 
given and the simultaneous equations are formed. Loading 
for these cases can be found by simply obtaining from charts 
the coefficients corresponding to the wing geometry, Mach 
number, and lift-curve slope, inserting in the appropriate 
equations and solving. 

Integration formulas for the loading distributions are 
given which enable the aerodynamic coefficients C t and C ti 
to be found. The rolling moment due to sideslip of a wdng 
with dihedral is shown to be equivalent to that of aileron 
deflection and a procedure for determining its value is given. 

For the special case of straight-tapered wings, the loading 
distributions and values of C t and Ci s are given in the chart 
form for a range of wing plan forms. 

Experimental and theoretical verification of the theory 
is shown to be good. The theory is applicable for large 
aerodynamic angles, provided the flow" remains unseparated. 
The compressibility considerations are reliable to the speed 
of sound subject to the limitations of the linearized com- 
pressible-flow equation. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., December 22 , 1949. 
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Fir, i re 15. Correlation of theoretical and low-speed experimental damping-in-roll coefficient Ci p for various plan forms 












APPENDIX A 

EQUATIONS FOR THE DETERMINATION OF ANTISYMMETRIC LOADING 


From NACA Rep. 921 (reference 2), the aerodynamic 


loading is obtained by solution of linear simultaneous 

equa- 

tions 


/ w\ m 

OC v ( y ) ^ j Avn G n } V = 1 , 2 

, . . . m 

(Al) 

where 



II 

ik 

» 


(A2) 

•A.n 2ft.. | ^ for ti — v 

,/b\ , . I 


(AS) 


UNSYMMETRIC LOADING 

by V and b vn are coefficients independent of plan form. 


gm ~ 


— 1 VLjy, 0 )/»» + -£ O'; M + l)./n.J/+l~t~ 


2(M+1) 


[ 


K=1 J 

(A 4 ) 


where 

J nfl are coefficients independent of plan form. 

The L (^m) functions, which can have ij v negative to find 
unsymmetrical loading, are given by 


for 0 or fx < 


M+ 1 


L(v, m)= Y 


v 1 + (o) tan +(-) 


+ 


and 

for < 0 or fx > 




,) t an A 


(c) {v ' 7? " ) 


M + 1 


H*, fO=i 


V 1 + (c ) tftn a (I»?»!+ 7 m)J +(jr) 

■(^“^ (I 1 ?* 1 ! 4 - v-) 


1 + 


tan A 


(e) (v ’~^ 


2 tan 

[i+w( 

B 

t an A 

‘<B 


[O 

(M- 

7j „) tan A 1 -f- 

(B 

1 (k»I+r 

h) ( a n aJ 


where 


VTT ... 

v*~ cos j spanwise position at which downwash is computed 


77^— cos -ng-pj spanwise position of incremental loading at the one-quarter-chord line. 


y (a 5) 


The above equations involve computations over the entire 
wing. However, if the loading is assumed to be symmetric 
or antisymmetric, the computations can be reduced to less 
than half the work. The case of symmetric loading is 
developed in reference 2 and the antisymmetric case is 
developed in the following section. 


ANTISYMMETRIC LOADING 

For antisymmetric loading, the loading on each side of the 
wing has the same magnitude and distribution but with 
opposite sign, or 


(X v OC, n f I — r 

G v — m + 1 — p j Ol G n = G m -j- 1 _ n 


(A 6) 


Equation (Al) can then be written as 

m — 1 
2 

oi v —^ (A vn — A Vtm ± i- n )G n (A 7 ) 

v= 1 

where the summation is only to — ? since ffmH=0 for 

* 2 

antisymmetric loading. 

With equations (A.‘ 3 ), equation (A 7 ) becomes 
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(The prime indicates the value for n = v is not summed.) 
Now, from equation (A4) 

— 1 [“ 0) (J no o) | 

&*n+g', »+*-- = 2r3/+i) L 2 


where 


and 


2 (M + 1 ) L 

L(v, A/jr 1) (J n. 3 fj-l — jm+1-n, -U-ul) 
2 

M "1 

L(v, n) (J nil Jm+l-r, 

M-i -I 

2 171 

/’«M= TT 2 Ml sin Mi0* COS M.tfv 

W m t mi = i 


+ 


(A9) 


(A 10) 


71 T J _ M *■ 

*"“m+T ^“M+l 

From equation (A9), /**“-./)« +i-«, m can be defined as 

,/ mi—jfui jf m f l-«, 

then, using equation (A 10), 

2 m , x 

/%„*= — S Ml COS ( Sin SIR 

i 1 ^i = l 
2 TO 

= - , , > ) Ml COS Mi 4 sin Mi</>«(f+COS Ml*) 

m + 1 “i 

and, since the terms of the summation for odd mi vanish, 

/* = — — X Mi sin AO0* cos mi</v (All) 

* m+ 1 Ml = 2,4,fi. . . even 

From* equation (All), 

J*nn— f*n,M+l-n (A12) 

Combining equation (A9) with (Al2) and defining 


ff v n (fvn [Jr, m+l — n 


then 


Af-f 1 

0 \ n = - - 1 -- [I(,.^) + L(,,i\/+l-M)]/%, (A13) 

2(M rl) m ~ o 

where for = 0 and ^ J* nv is equal to half the values 

given by equation (All) in order that the products can be 
fitted into the summation. With equation (A 13), equation 
(AS) can be written as 

m— 1 

a,=(2r ,+ -* ff *„) G'-IZ i 20 ™-}- »*«) (AU) 

-=1,2,3, . 


where 


m — 1 
* * ~2~ 

C „ = b yy b m _j_ I — P 

Ci >n~~~^vn by, — n 


From reference 2, 


b t 


sin <£„ 


n (cos <l> n — cos fa) 2 L 2 (m + 1) 


r i— ( — n-n 

2 L 2(m + l) J 


which gives zero values for for even (?? — v) values. Then, 
since m is odd, 

by, m -f- 1 — v ~ 0 


and 


C,= b w 


It should be noted that L(v.fx) simplifies somewhat for the 
antisymmetrically loaded wing since rj now is only positive 
in equation (A5). If only positive values of v are used, then 
equation (A5) can be written as 

L*(v, fi) = L* Vfl =L(rj , rj) + L(rj, — ij) 

= L(v , jx) + i(r, d/ -f 1 — m) 

In summary, the foregoing analysis for the antisymmetri- 
cally loaded wing gives 


m - 1 
2 

^ ) ftvnG n 
n=\ 


(A 15) 


where 


p=], 2, 3, . . . 


p vn ~2h vv -\-- g*y V for n=v 

c V 


= — 2 r,* + - - <7*„ for n 5* - 

C, 


Cv« byfj by m+l- n 

. OT + 1 


77? — 1 


(A 16) 


4 sin 


” (COS <t> 

g*„=!l*,n for n= v 


sin <ft„ T 1 —(—! )” ~ 

<j>„ — cos<£„) 2 |_ 2(m+l) _ 


Af+1 
' 2 


S *'* s= “2(M+l) S ^ 


/* — 


— XT S Ml Sin cos 
TYl T 1 mi « even 


T no -^forM = 0 


f * _T-foru — M+1 

/ „*±I— 2“ f01 2 
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For M=m, f* n „ simplifies to 

« _ p(— ])"+>■ sin 24> „~] / -sin 4<ft„ \ 

• "" L cos 20„ — cos 20 M \J -cos 40, 

fno—~ for n = 0 
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A=±i 




J* 


m + 1 




SVL i+ G0 < ’- ; ' ) "‘' ia J+(I) 


f 1 ) »7„) 2 1 j + 


Cv)(^ + ’ 7 ") 



i+l 

GO 

| (t? v — Vm) tan A 

+C0 

( '+ a (r) 

) r) v tan A 


0?» + 


vj 2 ) 

“T 


UmA^l+Q^tunAj+^y , 


1 +2 ^ ton A 


» 7 „=cos <£,, wlierc <#>»=- 


m + 1 

rj » — cos </> M where <£, — 


V n— cos <f>„ where <j> tl = 


M+ 1 

7? 7T 


m + 1 


For a discussion of the relative accuracies obtained for a 
choice of values of A/ and //?, see reference 7. The most 
favorable application is with M=m. 


With appendix A, the antisymmetrical loading on a plan 
form for any antisymmetrical distribution of a, can be found. 
The principal work in the computations is to obtain the co- 
efficients of the simultaneous equations (A15). These co- 
efficients can be presented in charts for the complete range 
of geometric plan-form parameters into which are introduced 
the effects of compressibility and section lift-curve slope. 
With the loading due to rolling known, the coefficients and 
derivatives are obtained by integration formulas. 

Section lift-curve-slope effect. — For a two-dimensional 
wing with the loaded line at the quarter-chord position, the 
position x aft of the loaded line where the induced downw ash 
equals the angle of attack of the wing can be obtained by the 
Biot Savart Law as 


CjcV 

2 



w= 

where T c 

or 


2 7T7 


w 

CiC 

or 

V~ 

4x7 a 



4x7 

then 

d a 

c 


7~ 

( c \ dci, 
\4x/ da 


APPENDIX B 

DERIVATION OF RELATIONS USED IN THE METHOD 
application of appendix a 

form geometry in the (bjc,)L,„ functions only, oi p,„ is a 
function of bjc , and sweep angle. As previously' shown, 
6/c„ is also a function of the spanwise variation of section 
lift-curve slope and is effectively equivalent to 6 /k,c„ where 
k, is the ratio of section lift-curve slope for a given Mach 
number at span station v to 2ir//3. The p,„ coefficients can 
be plotted against 6 /m\ with sweep angle as a parameter; 
however, 6 /k,c, will vary from zero to very large values for 
a range of plan-form geometry, and the plots become un- 
wieldy. For a range of aspect ratio, the values of &/k,c„ 
are a maximum for the zero tapered wings when tj,> 0.5 
(provided plan-form edges are not concave) and a maximum 
for the inverse-tapered wings for 7)„<0.5. The ratio of 
biKj, for r?,>0.5 for any plan form to those of the zero 
tapered wing or the ratio of 6 /k,c, for ij,< 0.5 for any plan 
form to those of the inverse-tapered wing gives a geo- 
metric parameter for any plan form that has maximum 
values that depend only on aspect ratio. 

The chord distribution for straight -tapered wings is 
given bv 

b_ = x l(l+X) 

c,. 


where dejda is the section lift-curve slope. Two dimen- 
sional section compressibility effects that do not follow the 
Brand tl-Glauert rule can be given consideration by taking 
the ratio of (dc t jda) tm possible at a Mach number to 2 t/ 0. 
Let k be the ratio of the section lift-curve slope at a given 
Mach number to 2? r//3 or (dejda)^ mpressibie= 2 tk/0, then 

a = k(c/2) 

Then the induced angle, k(c/ 2) aft of the loaded line, is 
equal to the angle of attack of the wing. For k= 1, this 
is at the three-quarter-chord line. For section lift -curve 
slope less than 27 t, k is less than one and the downvvasli is 
equal to the angle of attack at some point between the 
one-quarter- and three-quarter-chord line. 

To take into account the section lift-curve-slope variation 
in the present theory, the downwasli must be found at a 
distance a(c/2) aft of the loaded line. From the formulas 
of the summation in appendix A, bjc, should be taken as 
bjK,c,, where k, is the ratio of section lift-curve slope for a 
given Mach number at span station v, to 2 *-//?. 

Derivation of parameters for p, n . — The p,„ coefficients, as 
defined by equation (Alb) in appendix A, depend on plan- 
30 


Then, for X = 0, 


and for X=1.5 


= 2[1-|t?,|(1-X)] 


b 

Ac, 


(Bl) 


'2(1— |i?.|) 


b 

Ac. 


(B2) 


(B3) 


'2(2 + \ v ,\) 

The ratio of 5 /k,c, to equations (B2) and (B:3) gives, 
respectively, a geometric parameter as 

b 


-»?,)( * ) for 0.5 < v ,< O 
(bj Ac v )\=o \h v i 9 / | 


b j k,c , 

(bjAc.jy 


2(2 + 7?,) ^ J ^ for 0 <0.5 j 


Y (B4) 


Let H , be defined as two-fiftlis times the values of equa- 
tion (B4) (the fraction two-fifths is introduced to give H, 
the approximate values of to simplify plotting pro- 
cedures), then adding effects of compressibility (see Dis- 
cussion section) 

where 

4(1 — r} v ) 

5 


(h=- 


for 0.5 < T7„<1 


= 4(2±tp) f or o < 77 , < 0.5 
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For tapered wings, TJ V simplifies to 


jr 2(1— ^)(l+X )/gvl\ 

' 5[1 — 17,(1— X)j V K, / f<)1 05 ^ r l>< 1 

2(2 + „,)(l+\ )/,+ l\ 

25[1- 17 ,( 1 -X)]( «, ) fm ° 


(BO) 


Plots of ]t vn against TJ , in the range of //, = 0 to 4 will 
£i yo Pvn coefficients for wings of any elionl distribution for 
aspect ratios up to 10 or 12 . 

Linear asymptotes of p vn . — For large values of H n the 
p vn functions become linearly proportional f to TI V . Since 
'this linear characteristic appears at relatively low values 
of TJ y , the simple linear relation between p vn and TI V is 
quite usable. 

The L* vn function of appendix A is multiplied by hjc ¥ 
and the product is linearized. 

A (£)-A + (i^vl + v)*' m A ~ 1 

1 f A 

2n tall 2 "" 


-vL 




tan A 


. . .for 77<17 




. . 1 A 

Sm A ~2, tan 2- 


’-vL 


■V’+(jf I -0’' 
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y (B7) 


. . . for ?/>77 


-(cobv +tim A )© _ i + i sinA ~ 


2 ~ f an ^ ■ ■ • for 7] = r) 


cos A \e 


2 2 sm A _ _ 

1 _ ... for 17 = 0 

v v 


With the values of equation (B7) substituted into equa- 
tion (A 16) the values of p vn for arbitrary sweep angle are 
obtained. Thus, for m~l , the following equation (BS) 
gives values for p vn as 

;, 1 ,= VcosT + 1 - 02G tan a)//, + 7.968- 1.494 sin A + 

0.014 tan A + 0.082 (- - v 1 + 0 00 1 6 tan* A\_ 

2 V tan A ) 


o.ocsf ^X 1 + 0 0177 tnn? A \ 

\ t an A / 

( 1 — x 1 +0.1 7 1 7 tan 2 a\ 
\ t a iTa “ ) 


Pv2- 


0.034 

/ 0.851 
\ cos A 


tan A 

-2.901 tan 


A )//,- 


3.138+1.080 sin A- 
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->'1+0.0016 tan 2 
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— a 1 + 0.1 71 7 tan 2 1 
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lhi = ( ”, + 0 • 4 6 4 tail a)// 3 - 1.570-0.389 sin A- 
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\ / 1 — vl +0.1717 tan 2 A\ 

0.082 tan '2 + 0.231^ tanA / 

(\ - Vl +0.1 993 tan^Tv 
0082 (— ; 


„ _("°^2 8 _0.164 tan A )// s +3.417 + 0.083 sin A + 
p \ cos A / 

\ ( \ — an A \ , 

0.1 97 tan ‘-—0.082^ " tan A ) + 


/ 1 _ Vl +0.0SSG tan 2 

° 1G3 (v — — uin ) 

(\ - V l +0.1993 tan z_ A\ 

°-°K — i^a — ; 
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(B8) 


Linear spanwise distribution of (/cc),/(fcc) ar . With the 
condition that the product of section lift-curve slope and 
wing chord varies linearly spanwise, then 


2b 


kC- 


A k [1 (*rl k r)M 

and equation (3) becomes 
II v = d v (fiA k ) 


{ 1 — 7/ V [ 1 — {*t! k r) A] } 


1 ~K 


2 { 1 — r \ „ [ 1 — (kt/*h) A] } 


(B9) 


2 b 


and since kc is linear 
then 


r«(l +A) 
26 




krCr [1 +(W #c r) A] 
.1 


between 0 and 1 , k„ is a curve in the region between the linear 
function and a constant. 

Equation (B10) is given by figure 2 for m = 7^ but with 
the ordinate given by the parameter — 


where A t is the aspect ratio based on the wing chord equal 
to kc. Tn equation (B9), II, is reduced to terms of two param- 
eters Expressions of A. in terms of aspect ratio for 
straight-tapered wings and the distribution of section lift- 
curve slope can be found. 

For straight -tapered wings 


k>A)/l + A 

and equation (B9) becomes 

r m i j+(<c r /K B )x (P10) 

IT, - (i, |_ (lCs _4- Kr X)/l + X J 2 { 1 - r, , [ 1 - (Kr/Kfi)X] 1 


(kc) v 1 — 7 ? „ [ 1 — (k T I k r ) A] 

K v ~ K R “ “ ' ' 

c p 


1 — n.(l— X) 


(3.4/f(K fi + kj.X)/(1 +X)] 

and the abscissa by (kt/kk)X. 

For the case of linear distribution of (kc), and straight- 
tapered wings for which the chord and section lift-curve 
slope can be specified in three parameters, the loading and 
associated aerodynamic characteristics can be presented for 
a range of the parameters As, /3/l/[(K fi +KrX)/(l + X)], and 

(ktI K r) A. 

INTEGRATION OF ANTISYMMETRIC LOADING 

Rolling-moment coefficient and derivatives.— Rolling- 
moment coefficient is given by 


The distribution of k for straight -tapered wings is given by 

(Bll) 


j' G (v)v<h 


(B 1 2) 


where 


?j = cos<£> 

which, by an integration formula, 




m + 1 


m 

2 jivn) sin <p n 


(B 13) 


becomes 


^-0T) S «• “ s si ” *■ 


s I7^Tn S °n sin 2 4>n 

4 (m+1) Szi 


G 

Since the loading is antisymmetric, — | ±i=0 > aiwl 




(Bl4) 


For spanwise loading due to rolling, the loading is found 
as a function of pbj'lV, then equation (B14) divided bj 
pbj2V gives 




2 G n sin 2<f> n 


where 


lv 2 (m + 1) t = i 
G= Gf(pbj2V) 


(B15) 


The rolling moment due to ailerons will be found in 

appendix C. . 

Induced drag. The induced drag coefficient is, with equa- 
tion (B13), given by 


Equation (B10) is in terms of two parameters given by 
+ 1 I a nd (K r /K S )X. Solutions for spanwise load- 

L( k k+ K rX)/l +XJ . , , 

in* in terms of these two parameters and A* are valid toi Uu 
distribution of section lift -curve slope given by equation 
(Bl 1). Equation (Bll) indicates that at X= 1, k, is a linear 
function and at X=0, k, is a constant. For values of X 


pC D =pA £ G.a, w sin 4>, 


where a, is one-half the induced angle of the wing wake 
given by equation (Al4) for e,= », then for antisymmetric 
loading 


m — 1 


771 — ] 
2 


pC„ ( b„G*~ G, S'G.0«) s+ <t>, (Bl 6) 

H '* m + l7T\\ »=i / 
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For m = 7, tlie a Mi coefficients are equal to (for even mi) 


where the prime indicates the value of n=v is not summed. 

Span wise center of pressure, — The center of pressure on 
the wing half panel is given by 


Vc . p 


( l (1vdv 

_.'o 

f^Gdv 


The numerator is equal to $Ci/$A. If the Fourier series 
for loading is assumed, 

m 

<7(0) = S a M sin mi <t > , the denominator becomes 

Mi = even 

Mi/ Ml \ 

— I ) 2 \M! 2 — TM 


M i = 

then 


(B 1 7) 


f 2 . . ™ 

2_» <L, =sm 0 sin ]T) — a„,(— 1) 

i — pven */ o M,=fvon * 

- - 0O 

Vc . p . 

m Mi / \ 

Mt=fvcn \Mi — I / 

whore are the Fourier coefficients. 

Loading-due-to-rolling function and interpolation table. — 

Tlie Fourier series that approximates the antisymmetric 
loading with only a few terms is given by 


G(<t>) = a n Mi 4> 


(B18) 


The loading ft, is determined at span positions of 5j = cos 0„ 
= n^Fl ' T,IC arfi givon b -> r 


where <f>„ - 


2 C 1 

a H = ~ J o G{<t>) sin mi 4>d<t> (Bl9) 

With the quadrature formula of equation (Bl3), equation 
(Bl9) becomes, for antisymmetric loading, 


771 — 1 

4 “T" 


a "i = ^PY S ft. sin mi <t> n 


(B20) 


a 2 ; 


2 \ 2 


2 ft + ft+^j- ft 


■)1 


^4=2 (^i — ft) 


a *~\ 0? 6r ‘- <7 *+ > 2“ j 


(B21) 


Equation (Bl8) with (B21) can he arranged to give 


G {<!>)- 


1 ( \ 9 '9 \ 

= 2 ( sin 20 + sin 40+-^ sin (50 j (7,+ 

i (sin 20— sin C>0)ft + 

^ (t^- sin 20— sin 40 + ^ sin 60^ ft 


(B22) 


With equation (B22) the loading due to rolling can be deter- 
mined at any span position. Letting 0=0,.=^ and tabu- 

8 

lating the factors of G n as e nk} an interpolation table may he 
obtained to determine loading at span station k\ 

TABLE l\[,e nk 

Im- 7] 



0. 981 

0.831 

0. 550 

ft. 195 

/ 

/ 




7/2 

\ 


3/2 

5/2 

1 

0. 8155 

0. 5449 

-0. 1022 

0. 1084 

2 

— . 2706 

. am 

. 6533 

— . 2706 

3 .... 

. 1084 

-. ir.22 

.5449 

.8155 


ft — 2-j ?HkG n 

77= 1 


(B23) 


Equation (B23) may he used for interpolation of any form of 
loading coefficient, thus 

(B24) 

\CiC av / k 71 = I \CiC av / n 


V 



appendix c 

DETERMINATION OF ANTISYMMETRIC WING TWIST FOR FINDING SPANWISE LOADING DUE TO AILERON DEFLECTION 


WING TWIST FOR A GIVEN AILERON SPAN 

The determination of loading for an angle-of-attack 
distribution that contains a discontinuity by a method which 
satisfies the boundary conditions at a finite number of points 
can be made by increasing the number of points until the 
solutions become sufficiently accurate. For the method as 
given in appendix A, the number of points that satisfy the 
boundary conditions is given by m. For tbo huge \aluo o 
m required for accurate results, the computations become 
exceedingly laborious ; however, a procedure using a moderate 
value of m can be determined by use of a low-aspect -ratio 
theory with which a wing twist can be found that duplicates 
the results of the discontinuous angle-of-attack distribution. 

A theoretical but relatively simple method of finding 
spanwise loading due to inboard and outboard ailerons for 
wings of low aspect ratio is given by reference 4. In the 
present theory, as aspect ratio approaches zero, g* r „ values 
of appendix A become zero and the p,„ coefficients given by 
equation (AlO) become constant or independent of plan-form 
shape and equal to 

Vvn - n ) (C t ) 

p vv =2b v , ) 

These coefficients are given by the relations under equation 
(A15) and can be tabulated. 

TABLE Cl. pm 

[For 771 — 7 and **i=0] 


' Ss 'v v 

1 

2 

n 




1 

10. 1524 

-2.0000 

2 

-3. 0954 

5. 05C8 

3 

0 

j -2.0000 


0 

-1.5308 
I. 329<’> 


distribution given by reference 4 for a given aileron span. 

The aileron spans are arbitrarily chosen for the present 
theory as the mean value of the spanwise trigonometric 
coordinate of the down wash point at a section angle of 
attack equal to zero. For m=7, three aileron spans can be 
defined for both outboard and inboard ailerons. Let 17* be 
the aileron span, and 8 the spanwise point of the end of 
the aileron, then 

^=1 —cos 6 for outboard ailerons 
7j„=cos 6 for inboard ailerons 

For the present theory, the aileron spans defined are tabu- 
lated as follows: 

TABLE C2 


With equation (At 5), antisymmetric loading can be found 
for zero-aspect-ratio wings. As a comment on the accu- 
racy of the present theory for m = 7, the solution of equation 
(A 15), with the A = 0 p„ values for loading due to rolling 
crave the same values at the three semispan stations as does 
, (pbj2V) sin 2<t> 

reference 4, namely, — — 4 

The zero-aspect -ratio theory of reference 4 shows that 
all span loading characteristics are independent of plan- 
form shape for zero aspect ratio. This independence makes 
that theory ideal for obtaining the boundary conditions o 
the present theory for zero aspect ratio, which should apply 
with the present theory for higher aspect ratios for which 
plan-form shape has an effect on spanwise loading 1 he 
boundary conditions of the present theory are given by the 
antisymmetric values of <*„ in equation (A 1 5) . T ic pio ) etn 
is to find what antisymmetrical distribution of «> 1S rR_ 
quired for the present theory to duplicate the exact loading 
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Outboard 

Inboard 


Cast* 

T 

II 

TIT 

IV 

V 

VT 



57T 

7 IT 

5w 

3* 

0 

6 

10 

To 

16 

10 

1C 


V* 

0. 1085 

0.4444 

0.8019 

0. 5556 

0.8315 

1.0000 


For the aileron spans listed in table C2, the exact span 
loading distribution can be found from reference 4 Wit i 
the p., values listed in table Cl and the exact values of 
C,. G% and Gz from reference 4, equation (Alo) gives the 
twist required for the present theory to give the loading 
distribution for each ease listed in table C2 or 


— = 1 0.4524 


.0954 


/GA 

V 5 ) 
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~8 

« 3 
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“-Kt) UC 568 (t)'' 2 (t) 

1.530s(^)+«.329s($) _ 


(C2) 


The spanwise loading distribution from reference 4 for out- 
board ailerons is given by 


rwi = i 

l_ 5 _| outboard ^ 


sm 


(cos 0— cos 0) In 


(cos <£ + cos 0) In 


d+<t> 


sm 


-<t> 


cos 


2 

e+4>" 


cos- 


0— <t> 


2 J 


(C3) 



THEORETICAL ANTISYMMETRIC SPAN LOADING FOR WINGS OF ARBITRARY PLAN FORM AT SUBSONIC SPEEDS 


For tho full-wing-span aileron, or i?„= 1 

£ 


r<m~ 

=- cos <j> In 

[ 1 +sin 

L <5 

|l« = 1 7 T 

| COS <j> 

vd ailerons, 

with the same 

value of 

raw i 

jrmi 

p?(*n 

L G J inboard L ^ | *i„ = 1 

L 8 J„ 


outboard 


(C4) 


(C5) 


W i th equations (C3), (C4), and (C5), the spanwisc loading 
-> an(l at s P ftn stations <*>=tt/8, t/ 4, and 3 tt/ 8, or 
tj = 0.9239, 0.7071, and 0.3827 ran be tabulated for each of 
the cases given in table C2. 

TABLE C3 


o n 

6 

Case 1 

n 

III 

IV 

V 

VI 

Gi 

& 

0. 1136 

0. 1919 

0.2316 

0. 0451 

0. 1237 

0. 2373 

c 2 

0.500 






6 

. 2800 

.3851 

| . 1164 

. 3464 

. 3964 

u 3 







5 

.0190 

. 1022 

.3620 

.2922 

. 3754 

. 3944 


The twist distribution required for each case is obtained 
with equation (C2) and table C3, tabulating 

TABLE C4 


a 

7 

Case I 

IT 

TIT 

IV 

V 

VI 

on 

5 

1.0029 

0. 9713 

0. 9979 

0. 0444 

0.0128 

T. 0157 

<12 

5 

.0174 

.9957 

. 9913 

-.0169 

.9614 

.9788 

5 

.0056 

.0139 

. 9777 

1.0868 

1.0951 

1.1007 


— fo* * v iHjuuuon (A l o) 

ran be used to solve for spanwisc loading due to ailerons for 
any of the six eases. 


ROLLING MOMENT DUE TO AILERON DEFLECTION 

The rolling moment is given by 


A f*w 

Gi== ~4 J o G ( 0 ) sin 2 <p(l<p 


(C 6 ) 


For span loading due to ailerons, the loading distribution is 
distorted sufficiently such that the quadrature formula 
given by equation (Bl3) is not sufficiently accurate for 
m = 7 to integrate equation (C6). With equation (Bl8) 


r ttA 
e<=— a. 


(C7) 


Expanding equation (BlS) for tt/4, and 3W8, or 

obtaining G u G 2 , and G 3 in series of a’s, the sum of the G’s 
gives 


2 (°.7071<7,+£ i +0.707lG' 3 ) + a u — a lg -f-a 30 — a 3 


(C 8 ) 
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The higher harmonic coefficients can he put as factors of the 
G n . The rolling-moment coefficient becomes 


C r 


[ 


7 r 

16 


1.2071 G, 

(^14—^1 8 + ^30— a-34) 


0.707lTrr 

16 L + 


1.2071 G» " |(?2+ 

(^14— g 18 + fl3Q — a 64 ) 


l) ] 


’] 


G, 


(C9) 


1.207 lf? 3 

When h n is defined as the coefficients of G„ 

e« = A (hiGi + hiGj+k^G-i) (CIO) 

The ratio of («h— ffi 8 +a 30 — u- i4 ) to G n can be evaluated by 
the zero-aspect -rat 10 theory. It is expected this ratio will 
not vary appreciably with aspect ratio. From reference 4 
the loading series expansion gives for equation (Bl8) 

f<\ 4 

V 5 L, h j (cos 6 sin 11,9- 

/i, sin 6 cos (u,0) 


/ outboard TTjU, (/i, 2 — 

(<) 

\ 5 — 1 


+) =/M _(S) 

\ 8 /inboard \ 8 /, )o = , \ 5 /„ 


(Cll) 


^outboard 

These high harmonic coefficients are small, but are not 
negligible for loading due to ailerons. The h a are tabulated 
for each of (he rases 

TABLE Co 


fin 

C use I 

n 

III 

IA" 

V 

VI 

At 

hi 

hi 

0.1398 
. 1994 
. 1446 

0. 1388 
. 1963 
.1388 

0. 1379 
.1955 
. 1382 

1 0.1462 

.2004 
. 1400 

0. 1407 
.1973 
. 1391 

0.1402 
. 1975 
. 1397 


SPANWISE LOADING DISTRIBUTION 

The spanwisc loading distributions due to the twist 
distributions of table C4 art' found at three span stations, 
and, since these loadings are not completely defined by a few 
terms of the assumed loading series, the values of loading at 
other span stations cannot be found accurately by direct 
use of equation (B23) and table Bl. For zero-aspect -ratio 
wings, the spanwisc loading distribution due to aileron 
deflection is given to all span stations by equation (C3). 
The loading distribution for other than zero-aspect -ratio 
wings will fluctuate about the value given by equation (C3) 
in a manner similar to the manner that loading due to rolling 
varies about the function sin 2 <p of zero-aspect -ratio theory' 
Since the interpolation table of equation (B23) applies oniy 
to loadings that vary about the function sin 2<f> the loading 
due to aileron deflection can be divided by the ratio of equa°- 
Uon (C3) to sin 2<j> and the resulting loading will be approx- 
imat el v given by sin 2<p. 

The zero-aspect-ratio values of equation (C3) can be 
tabulated as ratios of 

sin 2<j> 


Define 


7? G(^ n )lB 
n sin 2 <»„ 


(Cl 2) 
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The zero-aspect-ral io values of R„ can he tabulated for each 
aileron-span case considered. 

TABLE C6— Rn 


The interpolation series of equation (B23) becomes 
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Outboard 


Inboard 


Case 
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IT 

III 

IV 

V 

VI 

\ 

\ n. 

0. 1685 

0. 4444 

0. 8049 

0. 5556 

0.8315 

1.0000 
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s s 

1 

2 ---- 

3 - — 

0. 1607 
.0500 
.0269 

0. 2714 
.2800 
. 1445 

0.3275 
. 3851 
.5119 

0.0642 

.1164 

.4132 

0. 1749 
1 . 3464 

j . 5309 

0. 3358 
.3964 
. 5578 


(Cl 3) 


where G=C and e Hk are giten 
0 


are 


hy table El. W i ill Tt k tab- 

kw 

ulated, values of loading at span stations r, l; =cos -j- 
obtained. 

TABLE C7 . — Rt 



\case 

I 

! 

TIT 

IV 

V 

VI 

0. 981 
.831 
.556 
. 195 

1/2 

3/2 

5/2 

7/2 

0. 1738 
. 1056 
.0341 
. 0233 

0. 2063 
,2787 
. 2250 
. 1212 

0.3162 
.3499 
. 43G5 
. 5304 

0. 0559 
.0802 
.2424 
.6363 

0. 1484 
.2533 
. 4225 

. 7343 

0. 3222 
.3589 
.4566 
. 7575 
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